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Transport in dimerized and frustrated spin systems 
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We analyze the Drude weight for both spin and thermal transport of one-dimensional spin-1/2 systems by 
means of exact diagonalization at finite temperatures. While the Drude weights are non-zero for finite systems, 
we find indications of a vanishing of the Drude weights in the thermodynamic limit for non-integrable models 
implying normal transport behavior. 
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Transport properties of one-dinicnsional spin- 
1/2 systems have recently attracted strong in- 
terest from the theoretical side (see e.g. |1I2I3I 
I4I5I6J and references therein). An intensively 
studied issue is the question under which con- 
ditions ballistic transport occurs at zero fre- 
quency and finite temperature T > 0, charac- 
terized by a non-zero Drude weight D. This 
quantity is the zero-frequency weight of the 
real part of the conductivity cr{uj), namely 
Rea{uj) = D{T)6{uj) -I- regular part. While it 
is known that the anisotropic Heisenberg chain 
has a finite Drude weight for thermal transport 
since the energy-current operator is a conserved 
quantity 0, there is a controversial discussion "5!, 
I3I4I5J whether transport remains ballistic if the 
model is extended by dimerization or frustration. 

In this paper, we will focus on two examples, 
namely thermal transport of dimerized chains and 
spin transport in frustrated chains. We consider 
the Hamiltonian (5 ~ 1/2) 
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H = J2hi = jJ2[XiSi- Si+i +aSi- S1+2]. (1) 
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Here, N is the number of sites and hi is the lo- 
cal energy density. We set A; = 1 if / is even 
and A; = A otherwise. Note that we use periodic 
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boundary conditions. 

The thermal Drude weight Dth{T) can be ob- 



tained from 7^ {pn — e 
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DthiT) = — 7-j Yl P" \{m\Jth\n)(" 



(2) 



while the Drude weight Ds(T) for spin transport 
follows fromHEl 



Ds{T) = 



N 



i-f) -2J2 pr 



l(™|js|n)|^ 



E„ 



E„ 



(3) 



jth and js denote the energy- and spin-current 
operator, respectively. They obey equations of 
continuity: i[H,hi[S[]] = -{jthls],i+i - jthls],i) 
with jth[s] = EiJth[s]j- The local energy den- 
sity hi is defined in Eq. Q and 5f is the local 
magnetization density. We refer the reader to 
Refs. 3 10 for full expressions of Jth[s]- Both 
current operators do not commute with H as 
soon as A 7^ 1 or a 7^ 0. The operator T — 

J J2i [Xi Si^Si'+i + 4a 5'i^5';+2 + H.c] is the kinetic 
energv P^llOj . 

Now we turn to the discussion of our numeri- 
cal results for the thermal Drude weight DthiT) 
of dimerized chains which we obtain by complete 
diagonalization of H. The data are shown in 
the main panel Fig. ^a) for A^ = 8, . . . , 16 sites. 
Since finite-size effects are strongest at low tem- 



peratures, we concentrate on T > J. Here, the 
leading contributions to Dth{T) take the form 



Ah(r)-ah,i/T^' + ah,2/r^ + 



(4) 



While we compute Cth.i directly by setting Z ~ 
2^ and e"''^" = 1 in Eq. 10), Cth,2 is extracted 
by a fit of Eq. Q to the numerical data at high 
T. As can be seen in the inset of Fig. ^a), both 
coefficients exhibit a significant decrease with sys- 
tem size though not monotonically. The system 
sizes are too small to draw conclusions about the 
actual finite-size dependence of -Dth, but the ob- 
served behavior does not support the notion of 
a finite Drude weight for TV — > oo contrasting a 
recent claim of Ref . . 

The second example is spin transport of frus- 
trated chains in the gapped regime. In the main 
panel of Fig. ^h), we show Ds{T) for a = 0.35 
and even A^ < 18. The main features are a fi- 
nite Drude weight at T = and a monotonic 
decrease of Ds{T) for T > 0.5 J. D^iT = 0) > 
has also been reported in Ref. JU] for a < 0.5 
and A^ < 20. Extension to larger systems at 
T — (e.g. with Lanczos techniques) should clar- 
ify whether Ds{T = 0) > is a property that 
survives in the thermodynamic limit A" — > oo. 

We further concentrate on the analysis at large 
T, namely the high-temperature prefactor Cg — 
lmix^oo[T ■ Dg{T)]. Its finite-size dependence is 
illustrated in the inset of Fig.^b) for a — 0.35, 1 
and A" = 8,9, ...,18. Apart from odd-even 
finite-size effects, which are more pronounced for 
a = 0.35, both data sets exhibit a monotonic de- 
crease of Cs with system size, as it is particu- 
larly obvious in the case of a = 1 . This provides 
strong evidence for Ds —> for A — > cxd in the 
gapped regime of frustrated chains which is con- 
sistent with Ref. 11 . 

In summary, we have numerically studied the 
thermal Drude weight of dimerized chains and the 
spin Drude weight of frustrated chains at T > 0. 
The finite-size analysis supports the conclusion of 
a zero Drude weight in the thermodynamic limit 
in these examples. We have found analogous re- 
sults for both kinds of transport for the models 
discussed here and spin ladders which will be re- 
ported in detail elsewhere (12,. 




Figure 1. (a) Main panel: Thermal Drude weight 
Dth{T) of a dimerized chain with A = 0.5 (see Eq. 
(JU). Inset: Cth,i,i = 1,2 vs. 1/A (see text), (b) 
Main panel: Drude weight Ds{T) of a frustrated 
chain with a = 0.35. Arrows indicate increasing 
A^. Inset: Cs vs. 1/A^ for A^ = 8, 9, . . . , 18. 
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